New Journal of Physics

Centrosome positioning is governed by mechanical interactions between the actin and microtubule cytoskeletons, yet the contribution of different actin architectures and their attachment to the cell membrane remains unclear. We use numerical simulations to study a centrosomal array of microtubules interacting with cortical, bulk, and centrosome-associated actin networks, including passive cross-linkers, motor proteins, and ERM-mediated coupling of actin to the cell membrane. We find that bulk and centrosome-associated actin alone are insufficient for robust centrosome centering, whereas sufficiently dense cortical actin promotes centering through steric interactions with growing microtubules. ERM-mediated attachment of cortical actin to the cell membrane enhances this effect, allowing stable centering at lower actin densities. In a two-centrosome system, microtubule interactions generate a spindle-like structure, while the organization of actin and its membrane attachment through ERM regulate centrosome separation. Antagonistic microtubule motors further limit centrosome separation and establish a finite spindle length. These results demonstrate that ERM-mediated coupling of actin to the cell membrane provides an important mechanical link between the actin cortex and microtubule-driven centrosome positioning, influencing both centrosome centering and spindle geometry.

Species at the highest trophic levels are generally the least resilient to perturbations, yet the factors that determine their persistence remain poorly understood. We ask how the intrinsic self-regulation of a basal host governs the persistence of the trophic levels above it. To address this question, we use a minimal discrete-time model of a host-parasitoid-hyperparasitoid system motivated by the biology of trigonalid wasps, obligate hyperparasitoids that complete development only when their caterpillar host is also attacked by a primary parasitoid. We compare two forms of host self-regulation in the absence of natural enemies: scramble (Ricker) fecundity, appropriate when the parasitoid is the host's principal enemy and which recovers the classical model of Beddington and Hammond as a special case; and contest (Beverton-Holt) fecundity, appropriate when the host is held in check by many additional enemies. Using numerical continuation, Lyapunov exponents, and direct computation of the model's trajectories, we find that the form of host regulation is decisive for persistence at the top trophic level. Under scramble host dynamics, coexistence is extremely fragile with respect to both parameter values and initial conditions. Under contest host dynamics, coexistence is robust across broad parameter ranges and loses stability only gradually, through a supercritical Neimark--Sacker bifurcation and a subsequent torus-doubling route to chaos. Therefore, in at least some trophic networks, the resilience of the highest levels can be determined at the base of the network.

We investigate the dynamics of the spread of the infection due to a phage- during the range expansion of an E coli colony. We find a rich panoply of spatio-temporal dynamics at large colony-level length scales, arising due to the interplay of mechanical forces, growth and viral infection dynamics at smaller cellular length scales. Strikingly, the phages "surf" the front of the growing bacterial colony by hitchhiking on E coli that are advected due to the growth of bacteria in the colony, resulting in an anisotropic spread of the phages. We identify microscopic processes - of the phage release during lysis (occurring on short millisecond timescales at length scales close to that of a single bacterium) and local nematic alignment of the rod-like E coli bacteria (occurring on timescales comparable to the bacterial growth rate and the length scales of a few bacterial cells) - that enhance the advective effects driving the hitchhiking behavior. Combining our experiments with a mathematical model, we explain the multiple spatio-temporal dynamical regimes - from the coexistence of uninfected bacteria, resistant (lysogenic) cells and phages, to fixed points where the entire population turns resistant or remains uninfected. Our work broaches a new frontier in the investigation of eco-evolutionary processes arising from the mechanical and ecological forces due to the interactions between viruses and their hosts.

Trigger waves are positive feedback-mediated traveling waves that propagate without slowing down or diminishing in amplitude, enabling biochemical coordination over distances where diffusion alone would be insufficient. Apoptosis can spread as trigger waves of caspase activity, as demonstrated in Xenopus eggs and egg extracts. We measured apoptotic wave speeds in extracts at two cell cycle states and found that waves propagated more slowly in M-phase arrested cytostatic factor (CSF) extracts than in interphase extracts (21 versus 27 m min-1). Luther's equation, v=2{surd}D{kappa}, predicts wave speed from the diffusion coefficient D and the positive-feedback doubling rate {kappa}. However, D and {kappa} were comparable between the two extracts and could not explain the speed difference. This discrepancy highlighted a key limitation of the classical equation: it does not account for an activation threshold in trigger wave initiation, and in fact the threshold required to trigger apoptosis was found to differ between cell cycle states. We therefore extended Luther's equation to include such a threshold, obtaining v{propto}{surd}(AD/Uc ), where A is the reaction rate and Uc is the activation threshold. Perturbing apoptotic sensitivity confirmed the predicted threshold-speed relationship. These findings establish the activation threshold as a key determinant of trigger wave speed and provide a broadly applicable scaling law for predicting the spread of switch-like biochemical signals.

Biological populations have been evolving in multi-dimensional trait spaces. Ecological interactions can induce the evolutionary diversification of a biological population into distinct multiple populations, which is called evolutionary branching in adaptive dynamics theory. If possible mutations are restricted to certain directions due to developmental, physiological, or physical constraints, the resulting evolutionary dynamics may be restricted to subspaces with fewer dimensions than the original trait spaces. In general, such constraints are also evolving slowly, due to rare and slight mutational deviations from the constraint subspaces. In this sense, constraints are not strict. However, the method for describing non-strict constraints that is globally consistent has yet to be developed. This study proposes a general and flexible method for constructing a non-strict constraint subspace from the corresponding strict constraint subspace. This method enables efficient analysis of adaptive evolution including evolutionary branching under non-strict constraints in the framework of adaptive dynamics theory. Applying this method to a resource competition model shows that the shape of the carrying capacity distribution can be treated as a non-strict constraint curve, and that the likelihood of evolutionary branching is affected not only by the shape of the constraint curve but also by its slow evolution.

Biflagellate zoospores are critical agents in the spread of Phytophthora spp., plant pathogens causing major agronomic and economic losses worldwide. These motile cells navigate soil water films toward plant roots using taxis triggered by environmental cues. Here, we investigate how potassium modulates zoospore motility by analyzing trajectory and flagellar beating parameters across KCl concentrations. Using high-speed imaging, we apply kymograph analysis and Fourier transform of the cell body direction over time to quantify flagellar dynamics, while low-speed imaging provides trajectory characteristics. Our results demonstrate that increasing potassium concentration significantly reduces swimming speed while enhancing directional changes. We establish a direct correlation between altered swimming speed and changes in flagellar beating parameters. Our findings indicate that potassium differentially affects anterior and posterior flagella beating. Furthermore, potassium exerts multiple effects, influencing beat amplitude, wave propagation speed and direction- change frequency. This reveals potassium as a key modulator of zoospore motility, triggering a shift from exploratory swimming to confined, reoriented movement. By uncovering the biomechanical responses of zoospores to potassium, this study advances our understanding of how surrounding signals regulate pathogen dispersal in the environment.

Most animals move intermittently. Of these, animals that fly and swim move in two stages: during active propulsion, when they beat their wings or swing their tails, and right after, when they glide or coast through the air or water. The kinematics of active movements are typically well studied, but the coasting period afterwards is less understood. We used larval zebrafish as a model of intermittent locomotion. Larval zebrafish swim in discrete bouts lasting ~200 ms, coast for short durations and distances, and then are largely stationary for ~800 ms. We find that larval zebrafish do not coast like a passive object. While they coast for distances that are correlated with the prior swim's strength, as expected from physics, they seem to be able to modulate the mechanics of their coasting: the time constant is not a constant, but varies from coast-to-coast. Furthermore, the time constant is correlated with the duration of the interval between swim bouts. This suggests that larval zebrafish can either shorten or extend their coasting duration, and that coasting inhibits the initiation of the next swim. Finally, the time constant changes quickly, while other swimming kinematics (e.g., distance swum) are autocorrelated over the course of at least several seconds. Altogether, this suggests that fish, and perhaps other animals that move intermittently, actively modulate their post-propulsion movement, possibly to gather sensory feedback in the absence of blurry visual feedback associated with fast locomotion.

We develop a general geometrical framework for growth processes in which the growth rate depends exponentially on the number of active growth-supporting components. The framework separates the fraction of active components from the size of the accessible component space, allowing qualitatively different growth laws to arise from the same underlying construction. Bacterial population growth provides an example of a fixed-size component space controlled primarily by activation and deactivation, whereas examples such as human population growth or speciation provide an expanding component space in which interactions between individuals promote further growth. This motivates a Replicative geometrical law, dX/dt= X exp({beta} X^{gamma}), whose exponent characterizes the scaling of the effective growth-supporting component space and is tested against historical human-population regimes and eukaryotic lineage-through-time data. For innovation indicators such as patents and arXiv manuscript submissions, where an existing item need not itself reproduce, the same construction motivates a Direct Geometrical production law without the prefactor X. The framework therefore provides a common description of fixed and expanding component spaces while distinguishing Replicative from Direct Geometrical growth.

Cells use ATP-driven chaperones to extract clients from biomolecular condensates selectively, yet no closed-form criterion links molecular parameters to condensate quality control. We derive a minimal non-equilibrium sorting reactor model that reduces the reaction-diffusion dynamics to one algebraic equation. In the idealized low-load limit, a concentration-dependent partition coefficient produces a saddle-node bifurcation between functional clearance and pathological storage. The critical Damk\"ohler number is \(\mathrm{Da}_{0c}=e^{-2}\) and is independent of the feedback strength. Corrections shift the value but preserve the topology. The model predicts a Kramers-type functional lifetime and recovery costs that scale as \(R^3\). Critical slowing down near the bifurcation provides early-warning signatures for FCS and FRAP.

Clathrin-mediated endocytosis (CME) is the primary pathway for internalization of extracellular and membrane cargo in eukaryotic cells. It is characterized by a patch of the plasma membrane invaginating and pinching off to generate a cytoplasmic vesicle. Recently, experimental evidence has begun to accumulate in support of proteins assembling into liquid-like condensates at nascent CME sites, but potential functional consequences are not understood. Given that condensate capillary forces can influence membrane bending, we constructed complementary analytical and numerical model systems to probe the mechanical influence of a protein droplet on endocytosis. Model evaluations and analysis reveal that a droplet can regulate the progression of endocytosis by stalling invagination until the droplet grows past a threshold size. After the droplet-mediated brake is released, the droplet provides a force that can drive the membrane bud toward scission. Thus, we propose a ``mechanical timer'' mechanism as a function for protein phase separation in cellular membrane trafficking.

Regular environmental cycles such as day and night or the lunar cycle are ubiquitous and drive many biological processes. Many organisms have evolved endogenous time-keeping mechanisms to anticipate these cycles to accurately time their behavior or synchronize physiological processes. The evolution of such synchronization to environmental cycles and their variability is not very well understood. As an example, we use the marine insect Clunio (Diptera, Chironomidae), which synchronize adult emergence and reproduction with the lunar cycle, and show different emergence patterns across populations. We use a model constructed with replicator-mutator dynamics to provide a novel outlook at the evolution of these emergence patterns. In particular, we investigate the effect of two traits of individual phenotypic variability of emergence, the spread and shape, on population-wide emergence patterns. We find that the semilunar and lunar emergence patterns are alternatively stable and occur when the phenotypic distribution has a narrow and is roughly flat. Arrhythmic emergence occurs when the phenotypic distribution is wide and nearly flat. The semi-lunar pattern can only be observed when the parameter of rhythmicity is changed. With our model, we introduce the concepts of frequency-dependent selection and individual phenotypic variability into the theory of chronobiology, which gives us insights into the evolution of rhythmicity in these organisms.

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

Adherent cells must continuously adapt to rapid environmental fluctuations to preserve mechanical integrity. However, conventional theoretical frameworks predominantly rely on quasi-static assumptions, thereby limiting their ability to capture the transient dynamics of cellular responses to high-rate perturbations, such as acute osmotic shocks. To address this limitation, we developed a biophysical model incorporating a Hill-type law that couples active cytoskeletal mechanics to cell-edge velocity, explaining the counterintuitive experimental observation that traction forces transiently decrease during rapid hypotonic swelling despite an increase in cell size. We further construct a kinematic-geometric phase diagram that connects classical quasi-static theories with our dynamic model. Moreover, we show that traction-force dynamics are sensitive to the loading rate and amplitude of osmotic shock. In addition, stiffness-associated simulations and supporting perturbation evidence suggest a cellular stiffness-associated recovery trend in which stiffer cells exhibit faster volume recovery driven by stronger recoil of hydrostatic pressure.

Vascular remodeling responds to local wall shear stress, but transport efficiency belongs to the network as a whole. We connect these scales by minimizing dissipation under a fixed maintenance budget, with vessel cost proportional to mass raised to an exponent b. At an optimum, dissipated power is proportional to vessel cost. This yields the shear set point |{tau} | {propto} r^(b-1){ell}^(b-1)/2 for radius r and length {ell}, with uniform shear recovered at b = 1. The associated local adaptation rule decreases the global functional, although the topology it reaches retains a memory of development. Simulations with 400 and 700 sinks show that staged growth improves transport mainly by allowing pruned connections to regrow. Fluctuating demand sustains loops without increasing the maintenance budget. The local condition fixes allocation on a given support, whereas growth and variable demand shape the topology that survives.

Phase separating biopolymers are often described by the sticker-spacer model, which accounts for non-attractive linkers that connect discrete attractive moieties. "Associative polymer" theories describe these systems as two coexisting fluids of stickers and spacers, respectively. In such theories the stickers and spacers are decoupled in the sense that sticker binding does not affect spacer conformations or vice versa. Informed by coarse-grained simulations, we develop a phenomenological mean-field theory to determine the microscopic coupling between the sticker and spacer degrees of freedom. Accounting for linkers connecting adjacent spacers introduces free energy contributions due to conformational entropy of the spacer and the binding combinatorics of the stickers. These contributions have opposite dependencies on spacer length because longer spacers allow a greater volume to search for binding partners but incur a greater entropy cost upon bond formation. The length dependent terms from these contributions cancel leaving a length-independent elastic term arising from the fact that a small amount of spacer stretching lowers the free energy by increasing the number of possible binding partners. The cancellation of spacer length from the theory justifies the decoupling approximation and supports mean-field treatments of heterogeneous spacers.

Behavioral variability reflects differences in internal and external parameters and provides a substrate for natural selection. However, uncovering its structure is challenging because behavior evolves nonlinearly across multiple timescales. Our goal is to reconstruct phenotypic variation from finite behavioral observations without prescribing the scales at which individuals should be compared. To achieve this, we introduce a multiscale dissimilarity that compares dynamics across all statistically resolvable resolutions by constructing reduced transfer operators at progressively finer scales and weighting their differences by finite-sampling uncertainty. This dissimilarity defines a geometry of dynamical phenotypes, which we reconstruct using diffusion maps. We show that the resulting phenotypic coordinates locally follow the parameter directions to which dynamics are most sensitive. We validate our approach in a stochastic model and apply it to bacterial and larval zebrafish behavior. In bacteria, we identify run speed as a major component of inter-individual variability, associated with the chemotaxis protein CheB. In zebrafish, we find that prior prey experience reshapes behavioral dynamics and phenotypic variability, with paramecia experience concentrating fish in a prey-capture-associated region of phenotypic space. Together, these results provide a framework for identifying internal and environmental variables that structure dynamical phenotypes across scales.

Astrocytes actively regulate synaptic transmission and neuronal excitability, yet their role in orchestrating macroscopic cortical network regimes and slow-wave oscillations remains an active area of reasearch. This study investigates how bidirectional neuron astrocyte interactions shape emergent population dynamics using a computational network model of excitatory and inhibitory neurons coupled to an astrocyte. The results identify astrocytic feedback topology, rather than astrocytic coupling strength alone, as a key determinant of emergent cortical network dynamics. By systematically dissecting pathway-specific connectivity, it has been shown that the neuronal population driving astrocytic activation and the neuronal population receiving gliotransmission jointly determine whether the network occupies asynchronous irregular (AI), synchronous irregular (SI), synchronous regular(SR), asynchronous regular(AR) or quiescent regimes.Directing gliotransmission selectively onto excitatory neurons consistently promotes population synchrony regardless of the population influencing astrocytic dynamics, whereas selective modulation of inhibitory interneurons induces network quiescence via strong suppression. Under dual-target gliotransmission, network synchrony is dictated by the population driving astrocytic dynamics: excitatory-only drive promotes synchrony, while combined or inhibitory-specific drive preserves asynchronous states. Furthermore, the model reveals that astrocytic signaling kinetics provide an additional temporal control mechanism that regulates the frequency and persistence of self sustained up states.

A major transition in the evolution of eukaryotic cells was phagocytosis: the ability of one cell to engulf another. In some cases, the engulfed cell persists, forming an endosymbiosis which is often mutually beneficial. However, it is unclear whether the interaction between the ancestors was also cooperative. There are clear benefits to developing the ability to eat a competitor, but can it ever be beneficial to eat a cooperator? We use a dynamical systems framework to identify the ecological conditions under which phagocytosis can evolve from a system of cooperating microbes.Beginning with a system of two cooperating species, we use invasion analysis to identify the parameter regimes under which a phagocytosing mutant can invade and the subsequent population dynamics. We identify a critical threshold that determines whether an invading predator will coexist with the cooperative community or replace members. Repeatedly invading with different predators, we find that the system saturates when there are three coexisting species, though this community does not have a single optimal composition. However, if we incorporate additional costs, we find there is a single evolutionarily stable community of two cooperators and a pure predator. Ultimately, we show that it is possible for phagocytosis to evolve from a cooperative community and, moreover, coexist with the original cooperating species.

Buruli ulcer (BU) is a neglected tropical skin disease caused by infection with Mycobacterium ulcerans. The mode of its transmission into the skin has not been fully clarified, but it is believed to involve mechanical transfer, perhaps by insect bites. M. ulcerans secretes a toxin called mycolactone which causes necrosis of the skin and formation of large skin lesions. Understanding the development of these skin lesions over time is crucial to understand the pathogenesis of BU. We formulate a mathematical model to describe the space-time dynamics of M. ulcerans bacteria, mycolactone and skin cell densities as well as fibrin deposition. We investigate the relative contributions to ulcer formation of direct cytotoxicity of mycolactone and ischaemia following the loss of coagulation control. Our results are consistent with existing biological evidence that higher mycolactone concentrations accelerate cell death. Furthermore, the model indicates that the rate of wound expansion is primarily controlled by the diffusion and net birth rate of bacteria, which sets the overall scaling for the speed of wound expansion, along with two dimensionless constants that define the shape of the trajectory of wound edge position against time, namely the ratio of the mycolactone diffusion coefficient to the bacterial diffusion coefficient, and the ratio of the mycolactone degradation rate to the bacterial net birth rate. When mycolactone diffuses significantly faster than bacteria spread or when mycolactone degrades more slowly than bacteria reproduce, the initial development of the wound is faster. In contrast, the long-term wound expansion speed appears to be relatively unchanged across a wide range of the dimensionless parameters.

Biofilm formation--the establishment of cellular communities within a self-secreted polymer matrix--is a behavior exhibited by almost all species of bacteria. Bacterial biofilms are well described as viscoelastic materials, with properties that grant physical advantages, such as elastic resilience and viscous adaptability. Additionally, many biofilm-forming species produce a subpopulation of highly resilient spores. In the model biofilm- and spore-forming species Bacillus subtilis, sporulation and matrix production are regulated by a common gene pathway, provoking the question of how sporulation influences biofilm viscoelastic properties. Investigating this point will facilitate the management of both salutary and deleterious biofilms, especially through manipulating their establishment and dispersion. Here, we investigate the interplay of sporulation and viscoelastic properties using rheological measurements of B. subtilis biofilms with varied matrix and spore production. We find that, to a significantly greater extent than matrix production, biofilm physical development is strongly altered by activation of sporulation through the spo0A regulation pathway. Spore-deficient strains initially establish physically robust biofilms with linear viscoelastic properties and recovery behavior distinct from spore-producing biofilms. However, spore-deficient biofilms fail to maintain their physical properties over time. In contrast, spore-producing wild-type and matrix-mutant biofilms exhibit more stable physical properties. Our results demonstrate that cell-level gene regulation is associated with macroscopic mechanical transitions in cellular communities, with sporulation granting not only a cellular-level survival advantage, but also greater community-level physical stability.